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Iron self-diffusion in nano-composite FeZr alloy has been investigated using neutron reflectom- 
etry technique as a function of applied compressive stress. A composite target of Fe+Zr and 
57 Fe+Zr was alternatively sputtered to deposit chemically homogeneous multilayer (CHM) struc- 
ture, [ natrura! Fe75Zr25/ 57 Fe75Zr 25 ]io. The multilayers were deposited on to a bent Si wafer using a 
3-point bending device. Post-deposition, the bending of the substrate was released which results in 
an applied compressive stress on to the multilayer. In the as-deposited state, the alloy multilayer 
forms an amorphous phase, which crystallizes into a nano-composite phase when heated at 373 K. 
Bragg peaks due to isotopic contrast were observed from CHM, when measured by neutron reflec- 
tivity, while x-ray reflectivity showed a pattern corresponding to a single layer. Self-diffusion of iron 
was measured with the decay of the intensities at the Bragg peaks in the neutron reflectivity pattern 
after thermal annealing at different temperatures. It was found that the self-diffusion of iron slows 
down with an increase in the strength of applied compressive stress. 

PACS numbers: 66.30.Fq 



I. INTRODUCTION 

During recent decades amorphous and nanocrystalline 
metals and alloys have been investigated as an im- 
portant class of materials with the possibility of tai- 
loring their properties over a wide range by control- 
ling particle size and morphologyiii2i&i More recently 
nano-composite alloys, in which nanocrystals are sur- 
rounded by an intergranular amorphous matrix, have at- 
tracted a great attention due to their interesting struc- 
tural^ and magnetic properties One of the favorable 
way to obtain a nano-composite alloy, is partial crys- 
tallization of the amorphous alloyi2i2ii£ The alloy struc- 
ture obtained above the primary crystallization temper- 
ature of the parent amorphous phase, but below the 
secondary crystallization temperature, has been termed 
as nano-composite phased Above the secondary crys- 
tallization temperature, the nano-composite structure 
fully crystallizes to form an equilibrium state of the al- 
loy. Structurally, the nanocrystals obtained after pri- 
mary crystallization are surrounded by an amorphous 
inter-granular phase to form a nano-composite phaseii^ 
Nano-composite alloys produced with an amorphous pre- 
cursor are the basis of interesting soft-magnetic alloys 
known as FINEMET 13 , NANOPERMii, HITPERMi^. 
Since the nano-composite phases produced in these al- 
loys is inherently a metastable phase, diffusion of the con- 
stituents would play an important role in understanding 
and determining their properties for long-standing appli- 
cations. Hence, atomic diffusion in such alloys is the key 
phenomenon for selecting their applications.- 1 - 

The situation becomes more complicated when the 
nano-composite alloys are produced in the form of a thin 
film. Deposition of thin films on to a substrate is known 
to produce films with a large intrinsic strain or stresses 
which often results from differences in thermal expansion 



(thermal stress) or from the microstructure of the de- 
posited film (intrinsic stress) ,iLi2iiS The intrinsic stresses 
may originate due to several factors (i) at the strained re- 
gions within the films e.g. grain-boundaries, dislocations, 
voids, impurities, etc. (ii) at the film/substrate inter- 
face due to lattice mismatch, different thermal expansion, 
etc. (iii) at the film/vacuum interfaces due to surface 
stress, adsorption, etc. or (iv) due to a dynamic pro- 
cesses e.g. re-crystallization, interdiffusion, etci^S These 
stresses may significantly affect the physical properties of 
the thin films, including atomic diffusion. 

It is known that when a material is deposited in the 
thin film state, the diffusion mechanism can be com- 
pletely different as compared to bulk state of that ma- 
terial, even when the material is in purely elemental 
form. Such a behavior has been mainly attributed to 
an increased defect concentration, metastability and un- 
relaxed state of the material. Therefore an extrapolation 
of bulk diffusivity may results in erroneous values of dif- 
fusivity in the case of thin films. Since many devices 
which are used for application are fabricated in the form 
of nm range thin films, self-diffusion measurements can 
be extremely important for their applications. 

In order to study the nature of stresses on self-diffusion 
we have chosen a simple binary FeZr alloy for this pur- 
pose. It was found that (as will be shown later) after an- 
nealing at 373 K, the alloy forms a nano-composite phase 
which further crystallizes above 600 K. The self-diffusion 
of iron was measured in the nano-composite state as a 
function of applied stress. The samples were deposited 
on to a substrate with a known bending. An external 
stress on to an [™ QtrW Fe75Zr25/ 57 Fe 7 5Zr 25 ]i multilayer 
was applied by releasing the bending of the substrate 
which resulted in an applied compressive stress on to 
the sample. Iron self-diffusion measurements were car- 
ried out using neutron reflectivity technique. It may be 
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noted that the neutron reflectivity is an excellent tech- 
nique for studying self-diffusion in nm range structures. 
Due to the fact that neutron scattering length densities 
for isotopes of an element are different, neutron reflec- 
tivity with depth resolution in sub nm range provides 
a unique opportunity for measuring self-diffusion. Con- 
ventional cross-sectioning and depth-profiling techniques, 
such as radiotracer, secondary ion mass spectroscopy 
(SIMS) are not suitable for measuring self-diffusion in nm 
range structures as the depth resolution available with 
cross-sectioning and depth-profiling techniques is of the 
order of a few nm. 

In an earlier worh^l we have demonstrated that neu- 
tron reflectivity is a technique which could be used to 
probe diffusion lengths of the order of 0.1 nm, and diffu- 
sion at temperatures less than 400 K could be measured. 
In the present work, the effect of compressive stress on 
the self-diffusion of iron in nano-composite multilayers 
was studied using neutron reflectivity. 



II. SAMPLE PREPARATION AND 
CHARACTERIZATION 



FeZr, CHM have been deposited on Si (100) substrates 
using magnetron sputtering technique. Small pieces of Zr 
rods were pasted on the natural ~Fe and 57 Fe targets in a 
symmetric way and the composite targets were sputtered 
alternately to prepare a chemically homogeneous struc- 
ture with nominal a composition Si/[™ at ™ raZ Fe75Zr25 (25 
nm)/ 57 Fe75Zr2s(10 nm)]io- The deposition of the multi- 
layer was carried out after obtaining a base pressure bet- 
ter than lxlO -6 mbar. During the deposition, pressure 
in the chamber was 5x 10~ 3 mbar due to 30 cm 3 /min Ar 
gas flow used for sputtering of the targets. All the sam- 
ples were deposited at a constant sputtering power of 50 
W. Before deposition the vacuum chamber was repeat- 
edly flushed with Ar gas so as to minimize the contam- 
ination of the remaining gases present in the chamber. 
Both the targets were pre-sputtered at least for 10 min. 
During the deposition the substrate was mounted on a 
specially designed 3-point Si wafer bending device. The 
substrate was oscillated with respect to central position 
of the target for better uniformity of the thickness of the 
deposited sample. 

In all the cases, thin Si wafers (300±10) /im were used 
as a substrate in order to avoid breaking during bend- 
ing. The Si wafer was fixed from both the ends, and by 
rotation of an asymmetric roller around the central axis, 
the bending height of the Si wafer can be varied between 
to 5 mm. A pin-lock system was incorporated so that 
release of bending by itself could be avoided. A com- 
pressive stress is applied on to the deposited film, when 
the bending of the Si wafer was released after deposition. 
The applied stress due to release of bending on to the 
Si wafer can be calculated using Stoney's formula^ and 
following a discussion given by Chen et &l~ The applied 
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FIG. 1: Schematic diagram of the bent Si wafer, used for 
calculation of radius of curvature. 



stress cr, is given by: 
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-, ) is the biaxial modulus of the silicon 

substrate and is equal to 180.5 GPa. Esi is Young's 
modulus for Si and vsi is Poisson's ratio for Si. Tsi is 
the thickness of the substrate, T/ is the thickness of the 
film and R is the radius of curvature. With the situation 
shown in the fig-H the radius of curvature can be written 
as: 
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Combining equation Q and J2J, the value of stress 
was calculated. The parameter used in the present case 
are T Si = (300±10)^m, T f = 370 nm, a = 40 mm and 
b was varied at 0, 3 and 5 mm. The obtained values 
of stress for the 3 cases are 0, 27 and 46 GPa. The er- 
rors in the calculation of applied stress were of the order 
of 15-20%, taking into account the uncertainties in the 
measured physical parameters. Samples with different 
known bending were deposited under similar deposition 
conditions. After deposition and release of bending, the 
surface profile of the samples was measured using a pro- 
filometer. It was found that the surface of the samples 
was flat and no changes in the surface profile were ob- 
served for a sample prepared with or without bending. 
This indicated that even after the bending the substrate 
gains its original state and the stress is applied on to the 
deposited multilayer. 

The composition of the deposited films was determined 
using x-ray photoelectron spectroscopy (XPS) depth pro- 
filing. The XPS profile was measured using monochro- 
matic Al K Q x-rays (1mm spot size) at the surface 
and at 3 different depths, after sputtering with Ar ions 
of 1 mA current and 3 kV accelerating voltage. The 
pressure during measurements was better than lxlO -8 
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mbar. The average composition of the films was equal to 
F e 75±3Zr25±3, excluding the data taken at the surface. 
Since at the surface, contributions from absorbed species 
like carbon and oxygen were significant, the average com- 
position of the film was determined with the data taken 
at 3 different depths. A small amount of oxygen was 
detected throughout the depth of the film. 

Structural characterizations of the samples were car- 
ried out with x-ray reflectivity (XRR) and grazing inci- 
dence diffraction using a standard x-ray diffractometer 
(XRD) with Cu-K a x-rays. The crystallization behavior 
of the multilayers was examined using differential scan- 
ning calorimetry (DSC) with NETZSCH, DSC equipped 
with extremely high sensitivity /i-Sensor. The conversion 
electron Mossbauer spectroscopy measurements (CEMS) 
were performed for determining the local environment of 
57 Fe in the samples. The measurements were carried out 
using a 50 mCi 57 Co-radiactive source in a Rh matrix 
and a gas flow proportional counter (He+4%CH4) for 
detection of conversion electrons. The isomer shifts were 
calibrated relative to a-Fe. Hysteresis loops as a function 
of azimuthal angle were measured using magneto optical 
Kerr effect (MOKE). 

Self-diffusion measurements were performed using neu- 
tron reflectometry technique at AMOR reflectometer at 
the Swiss spallation neutron source (SINQ), at Paul 
Scherrer Institute, Switzerland^ The reflectivity pat- 
tern was measured using two different angular settings 
(0.5° and 1.0°) in the time-of-flight mode. 



III. RESULTS AND DISCUSSION 



A. Structural Properties 
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FIG. 2: (color online) X-ray reflectivity of the as-deposited 
Si/[ natrm ' ai Fe 75 Zr25 (25±1 nm)/ 57 Fe7 5 Zr25(12±l nm)] 10 mul- 
tilayer at different applied stresses. The intensity shown on 
y-axis has been multiplied by a factor of 100, for clarity. 
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The multilayers prepared in this work have a period- 
icity only for iron isotopes, it is expected that x-ray re- 
flectivity of the multilayers would show a pattern corre- 
sponding to a single layer. Fig. [21 shows x-ray reflectivity 
pattern of the multilayer structure prepared at 0, 27 and 
46 GPa. As can be seen from the pattern, at the designed 
period of the multilayer there was no contrast for x-rays, 
which confirms the chemical homogeneity of the layers. 
The x-ray reflectivity pattern was fitted assuming a sin- 
gle layer and an 'oxide' layer of about 6 nm thickness on 
the surface of the multilayer, using a computer program 
based on Parratt's formalism^ 5 Such an oxide layer on 
the surface of the sample may stem because of absorbed 
oxygen or other light elements, when the samples are ex- 
posed to the atmosphere after deposition. Presence of 
such a layer was also evident from XPS measurements. 

On the other hand the neutron reflectivity pattern 
(fig- 13 1 showed well-pronounced Bragg peaks arising due 
to isotopic contrast between natural p e anc | 57 Fe. As can 
be seen from the figure, the sample prepared without any 
stress showed rather asymmetric Bragg peaks, while for 
the samples prepared with an applied stress, the peaks 
were more symmetric. Such an asymmetry of the Bragg 



FIG. 3: (color online) Neutron reflectivity of the as-deposited 
Si/[ natruro! Fe 7 sZr2 5 (25±1 nm)/ 57 Fe 75 Zr25(12±l nm)]io mul- 
tilayer prepared with applied stresses of 0, 27 and 46 GPa. 
The intensity shown on y-axis has been multiplied by a factor 
of 100, for clarity. 



peaks may arise due to incorporation of some free volume 
which may result in some internal strain or stresses dur- 
ing the growth of the film and might cause an asymmetry 
in the scattering length density. For the samples which 
were prepared in the bent state, the release of bending re- 
sults into an applied external stress on to the multilayer 
which eventually results in annihilation of free volume. 
A more detailed discussion related to this issue is given 
in the next sections. 

The neutron reflectivity pattern was fitted using a com- 
puter program based on Parratt's formalism^ and it 
was found that the pattern could not be fitted assum- 
ing sharp interfaces; instead a thin inter-layer of thick- 
ness (1±0.5) nm with the mean scattering length density 
of the two layers had to be introduced as inter-diffused 
layer. This means that already at room temperature 
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there is some amount of interdiffusion in the multilayer. 
The fitted parameters gives the structure of the multi- 
layers: Si/[ notrW Fe 7 5Zr 2 5 (25±1 nm)/ 57 Fe 75 Zr 2 5(12±l 
nm)]io, which is close to the nominal structure. 



B. Crystallization Behavior 

Prior to diffusion measurements thermal stability of 
the samples was studied with grazing incidence x-ray 
diffraction. All the samples were annealed together in 
a vacuum furnace in the temperature range of 373-573 K 
with 100 K step. In the as-prepared state all the sam- 
ples show a diffuse maxima cantered around 26 — 44.6°, 
with a width of about 4-5° (see Fig. 0J, which means 
that the samples are x-ray amorphous in the as-prepared 
state. The width of the diffuse maxima is comparable 
to the iron based amorphous alloys^ The average inter- 
atomic distance can be estimated using the relation, 26 
a = 1.23A/2 sint9, where 8 is taken to be the angle at 
the center of the diffuse maxima, and the factor 1.23 is a 
geometric factor which rationalizes the nearest neighbor 
distance with the spacing between 'pseudo-close packed 
planes'. As shown in the fig. r4J with an increase in 
the applied stress the position of the amorphous max- 
ima shifts towards higher angle side indicating a decrease 
in the average inter-atomic distance as shown in the in- 
set of the figure. Such a decrease in the average inter- 
atomic distance may be caused due to applied compres- 
sive stress. After annealing at 373 K, the broad hump be- 
comes narrow (width ~1°) and a nano-composite struc- 
ture is found, as shown in fig. [SJ The peak shape from 
the nano-composite structure could be fitted only by de- 
convoluting it into two lines, one corresponding to the 
parent amorphous phase and the second to a nanocrys- 
talline bcc-Fe phase. The area ratio of amorphous phase 
as determined from the fitting of XRD data was in the 
range of 15-25%. A slight decrease in the area ratio of the 
amorphous phase was observed. On further annealing at 
473 and 573 K, no significant changes in the XRD pattern 
of the samples were observed as shown in fig. [5] After 
annealing at different temperatures the position of Bragg 
peak shifts towards higher angle indicating a further de- 
crease in the interatomic spacing. Such a decrease in 
interatomic spacing is related with structural relaxation 
and is a consequence of annihilation of free volume. The 
grain size of the nanocrystals was about 10 nm, which 
increases marginally with an increase in the annealing 
temperature as shown in fig. |SJ 

Crystallization behavior of the samples was also stud- 
ied using conversion electron Mossbauer spectroscopy 
(CEMS). Samples prepared with and without applied 
stress were annealed at high temperatures. Fig. \7\ com- 
pares CEMS pattern of a sample prepared with and with- 
out applied stress before and after annealing at 473 K. 
The CEMS patterns were fitted assuming a sextate due to 
Fe nanocrystals and a doublet due to amorphous phase. 
Even in the as-deposited state a sextate with hyperfine 
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FIG. 4: (color online) X-ray diffraction pattern of the isotopic 
multilayers in the as-deposited state. The measurements were 
carried out in the grazing incidence geometry using Cu-Ka x- 
rays. The inset in the figure shows the change in inter-atomic 
distance as a function of annealing temperature. 
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FIG. 5: (color online) Grazing incidence x-ray diffraction 
pattern of Si/[ natrurai Fe 7 5Zr2 5 (25±1 nm)/ 57 Fe 7 5Zr 25 (12±l 
nm)]io multilayer prepared with an applied stress of 0, 27 
and 46 GPa after annealing at different temperatures. Open 
circles represents the measured data and the solid lines are fit 
to them. 
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FIG. 6: (color online) Average grain size as a function of 
annealing temperature shown representatively for a sample 
prepared with an applied stress of 0, 27 and 46 GPa. 

field of about 10 Tesla was found to be present. How- 
ever, the area ratio of this sextate is very small. This 
indicates that there is a small amount of magnetic al- 
loy while most of the alloy is not ferromagnetic. As the 
samples were annealed the contribution of this magnetic 
phase increases indicating an enhancement in the volume 
fraction of nanocrystalline Fe in agreement with the XRD 
results. It may be noted that the hyperfine field after an- 
nealing remains in the range of 18-23 Tesla while that of 
pure Fe is 33.3 Tesla. The reduced magnetic moment 
could result due to some thermal fluctuations. Table Q] 
compares the fitted values for the different cases as shown 
in fig. □ 



TABLE I: Fitted CEMS parameters for the sample prepared 
with and without an applied stress in the as-deposited state 
and after annealing at 473 K. 



Sample 


Sample 


Average 


Average 




Condition 


Hyperfine 


Quadrupole 






Field (T) 


(mm s _1 ) 


GPa 


As-deposited 


10.8±0.5 T 


0.40±0.01 


Gpa 


473 K, 1 hour 


18.4±1.3 T 


0.43±0.01 


27 GPa 


As-deposited 


10.2±0.9 T 


0.44±0.01 


27 Gpa 


473 K, 1 hour 


23.4±0.8 T 


0.70±0.01 


46 GPa 


As-deposited 


13.0±0.2 T 


0.44±0.01 


46 Gpa 


473 K, 1 horn- 


22.5±1.5 T 


0.42±0.01 



MOKE measurements were also performed to under- 
stand the changes in the magnetic properties of the sys- 
tem during amorphous to nanocomposite phases. The 
sample prepared without an applied stress exhibited no 
anisotropy as a function of azimuthal angle in the coer- 
civity. While the samples prepared with an applied stress 
clearly show uniaxial anisotropy, which persists even af- 
ter the nanocrystallization of the films (see fig. |SJ. This 
confirms that the bending stress induced in the films per- 
sists even after nanocrystallization. 

Formation of a nano-composite phase after primary 
crystallization of the amorphous phase is a general phc- 
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FIG. 7: (color online) Conversion electron Mossbauer 
spectroscopy (CEMS) pattern of Si/[ not ™ mi Fe 75 Zr25 (25±1 
nm)/ 57 Fe75Zr25(12±l nm)]io multilayer prepared with an ap- 
plied stress of 0, 27 and 46 GPa in the as-deposited state and 
after annealing at 473 K. 



nomenon in amorphous alloys. Often, it was observed 
that amorphous binary alloys crystallize in two steps. 
The primary crystallization reaction of most of amor- 
phous alloys leads to an evolution of nanocrystalline mi- 
crostructures whereas the phase formed after the second 
stage results in an intermetallic compound along with 
the nanocrystalline phase. The nominal reaction for such 
crystallization process had been given as: amorphous — > 
cv+amorphous — > a + (3; where a is the primary phase 
that precipitates out from the amorphous matrix and (3 
is an intermetallic compound In the present case 
crystallization of the amorphous phase can be regarded 
as the primary crystallization process and as evident from 
the x-ray data, the amorphous phase co-exists along with 
grains of cv-Fe. However, the primary crystallization tem- 
perature for the present case was found to be very low as 
compared with Fe67Zr33 amorphous alloyJi Since in the 
present case for Fe75Zr25 alloy, the Zr content is slightly 
lower, a decrease in crystallization temperature is not 
very surprising. In order to further confirm the structure 
of the alloy, a thin film with even lower Zr content was 
deposited under identical conditions of sputtering;— The 
composition of this film was FesnZr^o- The XRD pattern 
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FIG. 8: (Color online) MOKE measurements on Si/[' latrurai Fe 75 Zr2 5 (25±1 nm)/ 57 Fe 75 Zr25(12±l nm)]io multilayer prepared 
with an applied stress of 0, 27 and 46 GPa in the as-deposited state (a) and after annealing at 473 K (b). 



of this film showed a narrow peak even in the as- prepared 
state (not shown in the figure), indicating that the struc- 
ture forms a nanocrystalline state. Such a decrease in the 
primary crystallization temperature was also observed 
in an ion beam sputtered FessZris samplef^&Si and a 
phenomenon analogue to surface crystallization in amor- 
phous alloy ribbons22*2L22i2ii was found responsible for 
early crystallization of amorphous FessZris film in the 
thin film state. 

For the amorphous thin film formed by vapor depo- 
sition, the effective quenching rate is very high which 
results in a higher quenched-in free volume and results 
in an early crystallization of the amorphous phase as ob- 
served in the present case. In order to further understand 
the crystallization behavior of the alloy, DSC measure- 
ments were carried out under a constant heating rate of 
0.33 K/s. It was found that a very broad hump appears 
around 450 K and a relatively sharp peak appears around 
613 K in all the three samples. The hump appearing at 
450 K can be estimated as first crystallization step while 
relatively sharp peak corresponds to second crystalliza- 
tion step. From the XRD results, the first crystallization 
event was observed as early as 373 K, where the samples 
were annealed for 1 hour, in the DSC scan since the sam- 
ples were heated at a much faster rate, the crystalliza- 
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FIG. 9: (color online) DCS scans of the Si/[ natr " rai Fe 75 Zr25 
(25±1 nm)/ 57 Fe75Zr25(12±l nm)]io multilayer prepared with 
an applied stress of 0, 27 and 46 GPa. The shaded region 
shown in the figure correspondences to the temperature range 
used for diffusion measurements. 
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tion event was observed at higher temperature as shown 
in fig. [5J In a number of studies performed on bulk or 
thick films (thickness few jiva) , it has been found that the 
heat release during the first crystallization event is sig- 
nificantly smaller compared to the second crystallization 
event due to a slower diffusion at lower temperature^ In 
Al-based glasses, Foley et al2& have studied the crystal- 
lization behavior using transmission electron microscopy 
(TEM) and DSC. While growth of nanocrystals was con- 
firmed by TEM, there was no evidence of primary crys- 
tallization with DSC. Small diffusion of the constituents 
was argued for the observed behavior. In their case, they 
observed that, for diffusivity, D ~ lxlCT 19 m 2 s -1 , the 
level of heat output is nearly undetectable in DSC mea- 
surements, which requires a signal on the order of 0.1 
mW or greater. Unless the value of D was at least two 
orders of magnitude larger, the signal will be close to 
the noise level of the DSC. It may be noted that in the 
present case, the iron self-diffusivity around 400 K, is of 
the order of 10 -21 m 2 s _1 (as shown in later sections). 
Further, the grain-size as determined with XRD results 
was about 10 nm, the heat release for the formation of 
small grain sizes is expected to be small. In addition, the 
DSC measurements in the present case were performed 
in relatively thin film (370 nm) and the total mass ex- 
posed during DSC measurements was only 65 pgm which 
explains small heat release during the first crystalliza- 
tion event, in spite of the high sensitivity of the sensor 
used during DSC measurements. The presence of peak 
around 613 K, can be understood as second crystalliza- 
tion step. The onset of second crystallization tempera- 
ture was found around (608±5) K for all the 3 samples 
and there was no systematic effect of applied stress on 
the second crystallization temperature. 



C. Self-Diffusion Measurements - Time 
Dependence 

With the observed thermal behavior of the samples, 
for diffusion measurements, a temperature range for dif- 
fusion annealing was chosen from 413-533 K in order to 
study diffusion in the nano-composite state. This tem- 
perature region is also indicated in fig. [SI as shaded area. 
The three samples prepared with an applied compressive 
stress of 0, 27 and 46 GPa were first pre-annealed at 373 
K for 0.5 hour to obtain the nano-composite phase. For 
studying the time dependence of diffusivity, the samples 
were further annealed at 473 K and neutron reflectivity 
measurements were carried out after each annealing. In 
order to minimize the fluctuations due to a possible vari- 
ation in the temperature, all the samples were annealed 
simultaneously in the furnace. Fig. ^| shows a plot of 
neutron reflectivities as a function of annealing time at 
473 K. A relatively small time step was chosen in order 
to observe the structural relaxation of the samples. As 
can be seen from the figure, after annealing, the inten- 
sity at the Bragg peak decays. The initial decay was 



found to be much faster as compared to that with later 
annealing time. The decay of the Bragg peak intensity 
can be used to calculate the diffusion coefficient using the 
expression 36 : 

ln[I(t)/I(0)} = - 8ir 2 n 2 D(T)t/d 2 , (3) 

where 1(0) is the intensity before annealing and I(t) 
is the intensity after annealing time t at temperature T. 
The diffusion length Ld is related to diffusivity through 
the relation: 

L d = y/2D(T)t, (4) 

where t is the annealing time. The height of the Bragg 
peak was determined after subtracting the background 
due to Fresnel reflectivity by multiplying the data by a 
factor of q , where q is the momentum transfer. Fig. 1111 
shows an evolution of the diffusion length as a function 
of annealing time at 473 K. As can be seen from the 
figure, the diffusion lengths below an annealing time of 
600 s were found to increase much faster as compared 
to later annealing times. Such behavior in evolution of 
the diffusion length was also observed for Fe67Zr33 amor- 
phous sample 2 ™ and is a direct consequence of structural 
relaxation in the structuresi^S 

It is interesting to see that for the sample prepared 
without any stress, the diffusion length increased much 
faster as compared to the samples prepared with an ap- 
plied stress. As it is evident from the structural and 
magnetic measurements, samples prepared with an ap- 
plied stress exhibited a more relaxed state as compared 
to that obtained without an applied stress. It is expected 
that structural relaxation would be more dominant for 
the sample prepared without an applied stress. The over- 
all magnitude of the diffusion length follows the strength 
of applied stress and the degree of relaxation is propor- 
tional. 



D. Self-Diffusion Measurements - Temperature 
Dependence 

In order to measure the activation energy for diffusion, 
the samples were annealed in the temperature range of 
413-533 K with a step of 40 K. As can be seen from the 
El after an annealing time of 1800 s, in all cases, the 
fast relaxation process was almost completed, therefore 
for the calculation of the activation energy of the system 
all the samples were annealed for 1800 s at the above 
mentioned temperatures. It may be noted that annealing 
for 1800 s may not produce a fully relaxed state of the 
structure, even though for a comparison of diffusivitics 
for the samples prepared with different applied stress, 
the time for diffusion annealing was kept constant. 

Fig. ll2l shows the neutron reflectivity pattern obtained 
after annealing at different temperatures. Again in or- 
der to minimize a possible fluctuation in the temperature 
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FIG. 10: (color online) Decay of the Bragg peak intensity as 
a function of annealing time at 473 K for the sample prepared 
without an applied stress (a) with an applied stress of 27 GPa 
(b) and with an applied stress of 46 GPa (c). 

and annealing conditions, all the three samples were an- 
nealed simultaneously in the furnace. For a comparison, 
the neutron reflectivity pattern in the as-prepared state 
is also shown in the figure. As evident from XRD and 
CEMS measurements, after annealing the samples un- 
dergo from the amorphous to the nano-composite state, 
the intensity at the Bragg peaks increases marginally in 
the neutron reflectivity patterns. In a previous study, 
it was observed that at the event of primary crystalliza- 
tion in the amorphous Fe67Zr33 alloy 21 the neutron reflec- 
tivity pattern of an [ nat " ra *Fe 67 Zr 3 3(9 nm) 57 Fe 67 Zr 3 3(5 
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FIG. 11: (color online)Evolution of diffusion length 
as a function of annealing time and applied stress 
at 473 K in nanocomposite Si/[ natrura! Fe 7 5Zr25(25±l 
nm)/ 57 Fe 7 5Zr 25 (12±l nm)] 10 . 

nm)]2o multilayer, showed an increase in the intensity 
at the Bragg peak by a factor as high as 10, as com- 
pared to the as-deposited sample. Also the x-ray reflec- 
tivity pattern showed an appearance of a Bragg peak due 
to crystallization accompanied by a phase separation in 
the alloyiSi In the present case however, after anneal- 
ing at 373 K, the amorphous phase nanocrystallize, but 
the intensity at the Bragg peak increases only marginally 
(<10%). Also as shown in the fig. El 110 Bragg peak or 
structure due to a chemical period appeared up to 573 
K in the XRR pattern. This indicates that the primary 
crystallization behavior of Fe 7 sZr25 alloy is different as 
compared to that of previously studied Fe6 7 Zr33 alloy. 21 
However, since the matrix obtained after nanocrystalliza- 
tion showed no further significant changes between the 
temperature range 373-575 K (see also fig. (SJ, it is ex- 
pected that the diffusion process would be not interfered 
by structural changes. As shown in fig. 1121 the inten- 
sity at both the Bragg peaks decreases with increase in 
annealing temperature and after annealing at 533 K, the 
Bragg peak intensity almost vanishes. This indicates that 
after annealing at 533 K, both the natural and 57 Fe are 
layers almost completely diffused. With the procedure 
discussed in the previous section, the diffusivity at each 
temperature was obtained. Fig. 1141 shows a plot of diffu- 
sivities obtained with both the Bragg peaks for the three 
samples. As can be seen from the figure, both the Bragg 
peaks yield similar diffusivities within the experimental 
errors. The error bars in the present case are basically 
representing the errors in determining the height of the 
Bragg peaks obtained from a peak fitting procedure. 

The values for the diffusivity obtained for the three 
samples at the abovementioned temperatures (along with 
the separately annealed samples at 473 K), could be fitted 
to the relation, InD = IuDq — (E / 'k bT) , where Do, E and 
T are the pre-exponential factor, the activation energy 
and the annealing temperature respectively and is the 
Boltzmann constant. In all the three cases the logarithm 
of diffusivity follows Arrhenius behavior and accordingly, 
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FIG. 12: (color online) Decay ol the Bragg peak intensity as 
a function of annealing temperature for the sample prepared 
without an applied stress (a) with an applied stress of 27 GPa 
(b) and with an applied stress of 46 GPa (c). 



the activation energy and the pre-exponential factor for 
iron self-diffusion in Fe7sZr25 alloy was obtained. The ob- 
served values of both the E and Do are given in table [H] 
along with the values obtained for amorphous FegyZr33 
alloy. Fig. ^| shows a plot of diffusivities obtained from 
the 1st order Bragg peak for the sample at 0, 27 and 46 
GPa. The activation energy was found to increase with 
an increase in the strength of applied compressive stress 
(a steeper slope was observed with an increase in the 
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0.08 0.12 0.16 
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FIG. 13: (color online)X-ray reflectivity pattern of the 
Si/[ natr " ra! Fe r5 Zr 2 5 (25±1 nm)/ 57 Fe 75 Zr 2 5(12±l nm)] 10 mul- 
tilayers prepared with and without applied stress after anneal- 
ing at 373 and 473 K 



TABLE II: Activation energy and the pre-exponential factor 
for iron self-diffusion in nanocrystalline Si/[ natr " rai Fe 7 5Zr25 
(25±1 nm)/ 57 Fe75Zr25(12±l nm)]io multilayers as a function 
of applied stress during sample preparation. Both the activa- 
tion energy and the pre-factor represent the statistical aver- 
aged values obtained from the 1st and 2nd order Bragg peaks. 



Sample 


Sample 


Activation 


Prefactor 




condition 


energy 


(Do, mV 1 ) 






(E, eV) 




nano.-Fe75Zr25 


GPa 


0.24±0.05 


3xl(T 18±1 


nano.-Fe7 5 Zr25 


27 GPa 


0.31±0.05 


lxl0" 17±1 


nano.-Fe7 5 Zr25 


46 GPa 


0.51±0.05 


lxl(T 15±1 


amorphous- Fe67Zr33 


GPa 


0.38±0.05 


3xl0" 18±1 



applied stress). The result gives a clear indication that 
diffusivity for the sample prepared with applied stress is 
much slower as compared with that prepared without an 
applied stress. This result also supports the time depen- 
dence of the diffusivity as shown in fig. ^] 

In an earlier study Klugkist et al^L^i! studied Co and Zr 
self-diffusion in amorphous CoZr alloy using Radioactive 
tracer method as a function of pressure and temperature 
dependence. It was found that the pressure dependence 
for Co self-diffusion is extremely small while for Zr self- 
diffusion it is of the order of one activation volume. On 
the basis of obtained results it was concluded that Zr dif- 
fuses via thermal defects, whereas vacancy like thermal 
defects can be ruled out for Co self-diffusion. However, 
our results clearly indicates a decrease in diffusivity with 
an increase in applied stress. Here we would like to point 
out that a direct comparison between the studies per- 
formed earlier JSiii with that of our results could not be 
made. The following points are important to understand 
our results (i) The neutron refiectometry technique offers 
a depth resolution in sub nano meter range it is possible 
to measure initial stage of self-diffusion in a alloy unlike 
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FIG. 14: (color online) Arrhenius behavior of the diffusivity 
obtained from the 1st and 2nd order Bragg peaks for samples 
prepared with an applied stress of GPa (a), 27 GPa (b) and 
46 GPa (c). 



conventional techniques e.g. SIMS or radioactive tracer 
method (ii) As mentioned already that in our case, the 
alloy could have not attained a fully relaxes state as pre- 
annealing time was very short compared with that in 
literature. Combining the above mentioned points, it is 
not surprising that there is a strong dependence of self- 
diffusion of Fe on applied stress, which points out that in 
the initial state the diffusion mechanism could be differ- 
ent. 
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FIG. 15: (color online) Activation energy and pre-exponential 
factor for diffusion as a function of applied stress in 
Si/[ natrura! Fe 7 sZr 2 5 (25±1 nm)/ 57 Fe 75 Zr25(12±l nm)] 10 mul- 
tilayer. The data corresponds to the diffusivity obtained from 
the 1st order Bragg peak. The sample prepared with the 
highest applied gives and activation energy more than twice 
as compared with the sample prepared without an applied 
stress. The detailed values of activation energy and the pre- 
exponential factors are given in table ITTl 



Comparing the diffusivity for the sample prepared at 
GPa with that of amorphous Fe67Zr33 sample (also pre- 
pared at GPa), the diffusivity in the nano-composite 
state is slightly higher as compared to the amorphous 
sample. The activation energy for the nanocomposite 
sample was lower by 0.14 eV, while the pre-exponential 
factors were found to be exactly similar (see table ITTjl . 
An enhancement in diffusivity in the nano-composite 
state is not unexpected due to presence of grains and 
grain-boundaries (GB), while the amorphous phase is 
expected to be free from grains and GBs. However 
the enhancement in diffusivity in the present case is 
not as spectacular as observed e.g. in the FINEMET 
type nanocrystalline-Fe73.5Sii3.5BgNb3Cui in which the 
Fe self-diffusion showed a large enhancement over that 
in the parent amorphous phase4i*^ It may be noted 
that in the present case the composition of the nano- 
composite Fe75Zr25 alloy is also not similar to the amor- 
phous Fe67Zr33 alloy, therefore an enhancement in diffu- 
sivity may also occur due to increased concentration of 
Fe. In another study self-diffusion of iron was measured 
in the parent amorphous and nano-composite Fe^Zv^ 
alloy thin film produced by ion-beam sputtering. It was 
found that iron self-diffusion in both amorphous and 
nano-composite state was similar and found to occur ex- 
clusively through the GB regions which were amorphous 
in nature^ In the present case as well the GBs in the 
nano-composite state are amorphous, which might hap- 
pen due to the fact that in the nano-composite state, the 
structure consists of a mixture of nanocrystalline grains 
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of Fe and remaining amorphous phase. The nanocrys- 
talline grains of Fe would be surrounded by amorphous 
GBs and in such a situation a percolating path between 
the nanocrystals may not establish and diffusivity in the 
nano-composite phase would be similar to that in the 
amorphous state. 

On the other hand, the activation energy obtained 
for the sample prepared with highest applied stress was 
found to be larger (slower diffusivity) as compared to the 
amorphous Fe 6 7Zr 3 3 sample. This is somewhat surpris- 
ing as the sample prepared even at the highest stress is 
also in the nano-composite state. In case the diffusion 
mechanism is dominated by grains and GBs the effective 
applied stress should result in an enhancement of the 
diffusivity^ As discussed earlier, an applied compres- 
sive stress produced a more relaxed state of the sample 
as compared to samples prepared without stress. If the 
diffusion mechanism is dominated by a somewhat col- 
lective type migration of atoms both in amorphous and 
nano-composite case, annihilation of free volume would 
result into a diffusion mechanism involving a small group 
of atoms. 



IV. CONCLUSIONS 

In the present work, the effect of compressive stress on 
self-diffusion of iron in chemically homogenous multilay- 



ers of FeZr/ 57 FeZr was investigated. It was found that 
samples in the as-prepared state were amorphous and 
undergo primary crystallization when annealed at 373 
K. The diffusion measurements were performed in the 
nano-composite state and it was observed that with an 
increase in the strength of applied stress, the diffusivity 
decreases as compared to the sample prepared without 
an applied stress. An applied compressive stress on to 
the multilayer produced a more relaxed state of the sam- 
ple as seen from XRD. A diffusion mechanism involving 
a small group of atoms explains the observed diffusivity 
in the chemically homogeneous multilayers. 
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